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Abstract 

Using Finite-Size Scaling techniques we obtain accurate results for critical quanti- 
ties of the Ising model and the site percolation, in three dimensions. We pay special 
attention in parameterizing the corrections-to-scaling, what is necessary to put the 
systematic errors below the statistical ones. 

PACS: 75.40 Mg, 75.50.Lk. 05.50.-Hq, 75.40. Cx. 

1 Introduction 

The concept of Universality is perhaps one of the main discoveries of Modern Physics [Q . 
The critical exponents of phase transitions are among the most important quantities in 
Nature, as they offer the most direct test of Universality. Therefore, precise experimental 
measures of these exponents combined with accurate theoretical calculations are crucial 
cross-checks. Unfortunately, in three dimensions the range of variation of the exponents is 
very narrow. For instance, the correlation-length exponent, varies within a 10% interval 
for most systems |2|. Therefore, in order to distinguish between different universality classes, 
it is necessary to measure or calculate these quantities with several significant figures. 

There exist some powerful analytical techniques for computing critical exponents: e- 
expansions, high-T series, iV-expansions, or perturbative expansions at fixed dimension. A 
recent and complete study on this kind of calculations can be found in Ref. . A drawback 
of this approach is that the error estimate is quite involved. However, a 0.15% precision can 
be reached for v in Ising systems. 

A competing alternative is the use of Finite-Size Scaling (FSS) techniques Q combined 
with a Monte Carlo (MC) method, which in principle is able to measure with unlimited 
precision. The FSS method has the remarkable property of using the finite size effects 
to extract information about the critical properties of the system. In the language of the 
Renormalization Group (RG), we expect that for large enough lattices the divergences are 
fully described by the relevant operators. The MC method itself is not quite efficient as the 
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statistical errors in measures decrease only as the inverse square root of the numerical effort. 
However, the present sophisticated numerical techniques and algorithms, as well as the high 
computer power available, have allowed to reduce the statistical fluctuations largely. One 
could naively think that to get one more significant digit for a critical exponent is only a 
matter of multiplying by 100 the CPU time. This is not true in general, since the effects due 
to the finite-size of the simulated lattices eventually become larger than the statistical errors 
(in the RG language, the effects induced by the irrelevant couplings cannot be neglected 
anymore). Traditionally, one designs a simulation in order to get the systematic errors to 
lie below the statistical ones. With very high precision, a more quantitative treatment of 
systematic errors is required. 

In this work, we want to deal with the leading irrelevant terms (or corrections-to-scaling 
terms in the FSS language) in two of the simplest models in three dimensions: the Ising 
model and the site percolation 1^ . The reader might be surprised that the FSS ansatz holds 
for such a simple model as percolation, which is essentially not dynamical. The underlying 
reason is that bond-percolation is the g — > 1 limit of the g-states Potts model, as can be seen 
through the "Fortuin-Kasteleyn" representation of the latter The importance of both 
models has justified the construction of specific hardware as the Ising computer at Santa 
Barbara [Q, Percola § or the Cluster Processor However, the present update methods, 
as well as the power of the computers available allow to obtain very accurate measures for 
Ising models in general purpose computers. Regarding the percolation, an useful technical 
development has been the introduction of a reweighting method which allows to 

extrapolate the simulation results obtained at dilution p to a nearby p' dilution. As an 
outcome, dilution-derivatives can be also efficiently measured. This has suggested a different 
simulation strategy from the usual in percolation investigations ^, O, |l^, Instead of 
producing a small number of very large samples, we generate 0(10^ different samples in 
smaller lattices, in order to accurately measure derivatives with respect to the dilution and to 
obtain accurate extrapolations. The very nice agreement pc[ | with supposedly exact results 
for the critical exponents in two dimensions , and with other numerical results in three 
dimensions (see table ||), allows a great confidence in this new approach. In addition, the 
coincidence of two algorithmically different studies is a cross-check that reinforces both. 

The specific FSS method we use in this paper is based on comparison of measures taken 
in two different lattices at the value of the "temperature" for which the correlation-length in 
units of the lattice size is the same for both |l^ . Comparatively, this method is particu- 
larly well suited for the measure of magnetic critical exponents and for the parameterization 
of the effects induced by the irrelevant operators. We shall show that at the precision level 
we can reach (as small as 0.1% for the critical exponent i/, extracted from a given lattice 
pair), to take into account the effect of the leading irrelevant operator is unavoidable. For 
the two simple models we consider, very different situations are found. For site percolation, 
the scaling corrections exponent, w, is so large (w « 1.6) that other commonly ignored cor- 
rections, such as the induced by the analytic part of the free-energy, are of the same order. 
This makes our estimates of the critical exponents quite independent of the details of the 
infinite-volume extrapolation. But, on the other hand, the parameterization of the scaling 
corrections is remarkably difficult. On the contrary, for the Ising model we have oj « 0.8, 
and the infinite-volume extrapolation is mandatory. But the critical exponents related with 
higher-order corrections are large enough to allow for a neat, simple parameterization. 

In the next section we shall describe the FSS method we use. The measured observables 
are defined in section 3. The results for the Ising model and the site percolation are reported 
in sections 4 and 5, respectively. We will finish with the conclusions. 

2 Finite-Size Scaling 

Nowadays, a nice unifying picture of critical phenomena is provided by the Renormalization 
Group. In this frame, one can study not only the leading singularities defining the critical 
exponents, but also subdominant corrections (the Wegner confluent corrections In ad- 
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dition, from the Renormalization Group, a transparent derivation of the Finite-Size Scaling 
Ansatz (FSSA) follows (see Q and references therein). The starting point is the free energy 
of a d-dimensional system 

f{t,h,{u,}) ^ g{t,h,{u,}) + b-^U^^{byH,by-h,{u,byq), (i) 

where /sing is the so-called singular part, while g is an analytical function. We call b to the 
block size in the Renormalization Group Transformation (RGT), while yt, yh and yj {j > 3) 
are the eigenvalues of the RGT with scaling fields t, h and uj (j > 3). In the simplest 
applications (such as the ones we are considering) there are two relevant parameters: the 
"thermal field" , t, and the magnetic field, h {i.e. yt > Q,yh > 0) and we denote by {uj} the 
set of the irrelevant operators (0 > y^ > y^ > y^ > . . .) . One commonly uses the definitions 
v = 1/yi, T] = 2 + d—2yh and lu — —y^. The scaling field t can be identified with the reduced 
temperature in Ising systems, or with {p—Pc)/Pc in percolation problems. Taking derivatives 
of the free energy with respect to t or h it is possible to compute the critical behaviour of 
the different observables, including their scaling corrections [|l7| . A very similar strategy is 
followed in the study of a finite lattice, where we write for the free energy (see [|8| for a 
detailed presentation) 

fit, h, {u^}, L-i) = g{t, h, {«,}) + b-''f,,^g{by% by-h, {u^by^}, b/L). (2) 

At this point one takes b = L, thus arriving to a single-site lattice. By performing the 
appropriate derivatives, all the critical quantities can be computed. The result can be cast 
in general form for a quantity O diverging like t^'^" in the thermodynamical limit: 

0(i,t) =L^°/^ 

where Fq is a smooth scaling function. In usual applications one is interested in the ^ ^ L 
regime, thus is safely neglected. Of course in Eq. we have only kept the leading 
irrelevant eigenvalue, but, in fact, other scaling corrections like 

{Lyq,{Ly^+y^},... {t,j>3) (4) 

are to be expected. In addition, other kind of terms are induced by the analytical part of 
the free energy, g. For the susceptibility (or related quantities like the Binder cumulant or 
the correlation-length, see below) one should take the second derivative with respect to the 
magnetic field, h, in Eq. (|^). The leading contribution of the analytical part is independent 
of the lattice size, thus if one wants to cast the result as in Eq. (^), corrections like L^'*/'^ 
should be added. 

Equation (^ is still not convenient for a numerical study, because it contains not directly 
measurable quantities like ^{oo,t). Fortunately, it can be turned into an useful expression 
if a reasonable definition of the correlation length in a finite lattice, £^{L, t), is available: 

0{L,t) = L^°'' 

where Fq is a smooth function related with Fq and f^. 

To reduce the effect of the corrections-to-scaling terms, one could take measures only 
in large enough lattices. Even in the simplest models, as the two considered in this paper, 
if one wants to obtain very precise results, the lattice sizes required can be unreachable. 
However, we shall show that this is not the most efficient option. In the specific method we 
use, the scaling function is eliminated by taking measures of a given observable at the same 
temperature in two different lattice sizes (Li, L2)- At the temperature where the correlation 
lengths are in the ratio Li : L2, from Eq. (||) we can write the quotient of the measures of 
an observable, O, in both lattices as 

/L \^°/^ 

Qo\Q^^Li. = [j-^j + Aq^l:,'^ + . . . , (6) 



Fn 



L 



e(oo,t) 



+ o(L-",r'") 



(3) 



Fc 



(5) 
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where is a constant. 

The great advantage of Eq. (||) is that to obtain the temperature where — L1/L2, 
only two lattices are required, and a very accurate and statistically clean measure of that 
temperature can be taken. In addition, the statistical correlation between Qo and 
reduces the fluctuations. Other methods, such as measuring at the peak of some observable 
suffer in general from larger corrections-to-scaling. Computing the infinite volume critical 
temperature and measuring at that point performs well for studying observables that vary 
slowly at the critical point, as those used for computing the v exponent. However, the 
magnetic exponents require measuring quantities that change rapidly with the temperature 
and this is more involved. We think that our method outperforms any other previously 
used, specially in the computation of the rj exponent. 

To perform an extrapolation following Eq. (^, an estimate of ui is required. This can 
be obtained from the behaviour of dimensionless quantities, like the Binder cumulant or the 
correlation length in units of the lattice size, ^(i, t)/L, which remain bounded at the critical 
point although their i-derivatives diverge. For a generic dimensionless quantity, g, we shall 
have a crossing 

5(L,t"-°-(L,s)) = g(sL,t"°-(L,s)). 
The distance from the critical point, t'^''°^^(L,s), goes to zero as | |T9| |: 

<"°^^(L,s)cxiy^L— 1/- . (7) 

From Eq. (Q), a clean estimate of uj can be obtained provided that \yi \ — uj and ^ jv — lo are 
large enough (say of order one). 



3 The Models 

We will consider a cubic lattice with periodic boundary conditions and linear size L, the 
volume being V — . In the case of the Ising model we consider the usual Hamiltonian 

<ij> 

where the sum is extended over nearest neighbour sites and the spin variables are ±1. 
The fundamental observables we measure are the energy, and the magnetization 

<i,j> i 

The energy is extensively used for f3 extrapolation ||2^ and for calculating /3-derivatives 
through its connected correlation. 

The other quantities that we measure are related with the magnetization. In practice we 
are interested in mean values of even powers of the magnetization as the susceptibility 

X - V{M') , (10) 

or the Binder parameter 



54 - o - o /A.2\2 • ill) 



3 _ 1 (M^) 
2 2 {M^ 

The cumulant 54 tends to a finite and universal value at the critical point. As correlation- 
length in a finite lattice, we use a quantity that only involves second powers of the magne- 
tization, but uses the Fourier transform of the spin field 

5(fc) = -i^e"= V. ■ (12) 
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Defining 

= -^(|CT(27r/L, 0,0)1^ + permutations) , (13) 
we will use as correlation length 

The site percolation is defined by filling the nodes of a lattice with probability p. Once 
the lattice sites are filled (we call this particular choice a sample) a system of spins is placed 
in the occupied nodes. The spins interact with the Hamiltonian (^ at zero temperature 
(/3 = oo). In this way neighbouring spins should have the same sign, while the signs of spins 
belonging to different clusters {i.e. not connected through an occupied lattice path) are 
statistically uncorrelated. Thus, by counting the number of spins contained in each cluster, 
{ric}, we know the exact values of {M"^) and {M'^) in a particular sample: 

{M') = 4E-e, (15) 



V 

C. 



where the sums are extended to all the clusters. 

To compute the quantities involving Fourier transforms of the magnetization we measure 



^ic(fc) = -^Ee''= V. , (16) 

where the sum is extended to the sites of the c-th cluster, arriving to 

(|a(fc)p)=5]|n,(fc)p. (17) 

c 

We then average Eqs. (|l6| ) and ( |l7| ) in the different samples generated. This new average 
will be denoted by an overline. So we define the correlation length and the cumulant 174 as 



4sin^(7r/L) 

54 = 7: - o==2 • (19) 



2 2 



Another universal quantity, whose non-vanishing value proofs that the susceptibility is not 
a self-averaging quantity, is the cumulant (72 



^ ^-M^l^S^ . (20) 

A last technical comment for our percolation study is that we store the actual density 
values obtained with probability p, in order to perform a p-cxtrapolation of the mean values 
of the interesting observables, and also p-derivatives |lO|, |ll| , |2^ . 

Both for the Ising model and site percolation, the observables we use to compute the 
two independent critical exponents, 77 and z/, are 

X ^ x^y[2-ri), (21) 
dpS,, dp£, X = jy + 1 . 



5 



For the sake of completeness we will link here our method with the more classic approach 
followed in percolation. The basic entity in percolation is the number of clusters of size s 
divided by the lattice volume, Ug Q. This object induces a probability of finding a cluster 
of size s, given by sus- Near the percolation threshold, Pc, Us follows the law 

ns = s--f{s'/^p~Pc)), (22) 

where cr and r are critical exponents, and / is a scaling function. This yields just at p ~ pc 

ns = s-^ [A + Bs-^ + ...), (23) 

where J7 is a corrections-to-scaling exponent. We can relate the thermodynamical critical 
exponents, 77, v and w with the more standard exponents in percolation cr, r and fJ: 

^ = (24) 



ad ' 
(2-d)T + 3d-2 

av ' 



(25) 
(26) 



where d is the spatial dimension of the lattice (in opposition to the fractal dimension df 
(a^)-i). 



4 Results for the Ising model 

We have used a Single Cluster (SC) update algorithm ||2^ which is known to perform very 
well for this model. We take measures every 50 SC. We have accumulated 8 millions of 
measures for lattice sizes i=8, 12, 16, 24, 32, 48, 64, 96 and 128 at /3 = 0.22165. For 
the statistical analysis, we use a jack-knife method with 50 bins of data. Our pseudo 
random number generator has been a corrected shift register generator introduced in Ref. j24| 
improved by adding (modulus 1) a congruential generator (see Ref. [|o)). 

The results for the critical exponents using the quotients for the observables of Eq. ( |2l| ) 
from lattice pairs of sizes (L, 2L) measured where = 2 are shown in table |^. We also 
report the values of the universal 174 cumulant at the same points. From this table, it is 
apparent that, with the statistical error reached, an infinite volume extrapolation is needed. 
This is especially clear for the i] exponent. 



L 


h' 


77 


.94 


8 


0.64379(37) 


0.01097(40) 


0.72177(19) 


12 


0.63778(46) 


0.02094(31) 


0.71460(19) 


16 


0.63654(40) 


0.02548(38) 


0.71048(18) 


24 


0.63385(44) 


0.02927(34) 


0.70668(21) 


32 


0.63277(48) 


0.03129(38) 


0.70481(25) 


48 


0.63164(48) 


0.03273(37) 


0.70317(25) 


64 


0.6316(6) 


0.03376(39) 


0.70204(33) 


00 


0.6294(5)(5) 


0.0374(6)(6) 


0.6984(5)(6) 



Table 1: Critical exponents v and 77 for the Ising model obtained from pairs of lattices of sizes 
(i, 2L) where = 2. We present also g4{L) at the same points. The last row corresponds 
to an infinite volume extrapolation considering the leading scaling corrections. The second 
error bar is induced by the uncertainty in to: when uj increases, v and (74 increase, while 77 
decreases. 
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L 








8 




0.2216246(32) 


0.2218571(29) 


12 




0.2216438(13) 


0.2217279(18) 


16 




0.2216487(9) 


0.2216878(11) 


24 




0.2216519(5) 


0.2216661(7) 


32 




0.22165314(45) 


0.2216601(6) 


48 




0.22165438(23) 


0.22165709(30) 


64 




0.22165432(19) 


0.22165563(26) 


Table 2: 


Crossing points of ^/L and 54 for the Ising 


model obtained from lattice pairs 


(i,2i). 










xVd.o.f. 




/3c 


8 


10.50/10 


0.934(14) 


0.221654433(83)(18) 


12 


10.06/8 


0.938(24) 


0.221654447(90)(7) 


16 


1.796/6 


0.87(4) 


0.22165456fll)fl) 


24 


1.727/4 


0.86(9) 


0.22165459(15) (5) 



Table 3: Results of the infinite volume extrapolation of the crossing points for 34 and f /i, 
including data from L > Lmin in the Ising model. The first error bar in (3^ is statistical while 
the second is due to the uncertainty in uj. We quote our preferred value and error bars by 
underlines. The extrapolation for /3c decreases when uj increases. 



To perform this extrapolation, one should try to take into account the corrections-to- 
scaling. For the Ising model, we shall show that the leading order corrections are enough to 
obtain a fair extrapolation. 

We first need to evaluate the ui exponent. From the data of table it is difficult to obtain a 
sensible value of cu. We use the shift from the infinite volume critical coupling of the crossing 
points of the scaling quantities ^/L and 54 which are much more accurate. These points are 
shown in table |^. Then we carry out a joint fit of all values to the functional form given in 
Eq. (0) (see Refs. [|6[ |2|] for a more detailed exposition of the method). Our criterium is to 
fit the data, using the full covariance matrix, for lattices greater or equal than a given Lmin, 
increasing this minimum size until a stable value and a reasonable x^/d.o.f. are found (see 
table ^). We take as fitted parameters the corresponding to the first satisfactory imin with 
the statistical error of the fit discarding L^nin- In this case, we observe that Lmin — 16 is 
enough for our precision. The value obtained uj — 0.87(9) is compatible with the computed 
using analytical techniques |Q, or the recent experimental value uj = 0.91(14) Thus 
one can be confident that the correction-to-scaling are mainly due to the leading term when 
L > 16. We should remark that for obtaining uj an estimation of v must be used. For our 
accuracy, a value of i' with an error at the 1% level is enough. Therefore, we do not need a 
previous infinite volume extrapolation and it is safe to take v = 0.63 for this purpose. The 
obtained infinite volume critical point, /3c, is scarcely affected by the uncertainty in lo. 

Using this estimate for w we can perform a fit to the Eq. (H). In the last row of table 
we present the extrapolation results for i/, rj and (74. In all cases we have used the criterium 
described above, in order to deal with higher order scahng corrections (Lmin = 12, 16, 16 for 
ly, rj and 34, respectively). 

In Fig. |l| we show graphically the fit quality for the quotients used to obtain the exponents 
as well as for 174 cumulant. 

We remark that the pair (L, 2L) having a systematic error in rj smaller than the final 
error in our extrapolation has 2L = 2000 (2L = 800 for v). We recall that our largest pair 
has 2L = 128. 
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5.95 




Figure 1: Infinite volume extrapolation for Qogj 
lines correspond to fits with lo = 0.87 from L > L^m 



2^+-, Qap = 2-1-'', and 54- The solid 
. A dashed line is plotted for L < L^i^- 



Values for the critical quantities obtained with MC by other authors are reported in 
table ^. For comparison, a recent series computation[p6| yielded /3c = O.22I659I5. Only 
in Ref. an infinite-volume extrapolation is considered for the critical exponents. A 
different approach has been used in Ref. pof , where the Hamiltonian is numerically tuned, 
in order to make the 0{L~^) corrections for cumulant 174 vanishing. This largely reduces the 
corrections-to-scaling for the exponents. However, the data of have been analyzed as if 
these coefficients, Aq, would he exactly zero. However, there still is an error associated to 
the uncertainty in the assumption Aq = 0, that has not been considered in Ref. [ p9[ . In fact, 
there is nothing special in the value Aq — 0, the only essential ingredient for the dramatic 
reducing of the error estimate is to neglect the error in Aq ■ Had we disregarded the error 
in Aqq, we would obtain i' = 0.6294(2), 77 = 0.0374(2). However, we do not believe this to 



Ref. §c 1^ V d 7 

|2| 0.2216544(10) 

|18| 0.2216546(10) 0.6301(8) 0.037(3) 

|2| 0.221655(1)(1) 0.625(1) 0.025(6) 

ll" 0.2216544(6) 

0.6299(3) 0.0359(10) 

This work 0.22165456(15)(5) 0.6294(5)(5) 0.0374(6)(6) 0.3265(3)(1) 1.2353(11)(14 ) 

Table 4: Previous Monte Carlo determination of critical quantities for the Ising model. 
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be a valid procedure. 



5 Site percolation results 

The MC simulation in this case is rather different, since one generates directly independent 
configurations. We will work in the so-called canonical formulation in which the probability 
of finding a hole in a given lattice site is independent from the rest of sites. 

It is very fast to generate the different configurations and most of computer time is 
employed in tracing the clusters. We generate 32 millions of samples for L < 96, 16 millions 
for L = 96, 128 and 4 millions for L = 192. As we need individual measures for the 
^-extrapolation it is necessary to store them on disk as they are obtained from different 
processors. In all cases we simulate at p = 0.3116 [||. 

In table ^ we present the results for the exponents and r/ as well as the 52 and 54 
cumulants, obtained from different pairs of lattices. 

Regarding pc, we show in table || the crossing points of ^/L and (74 for pairs (L, 2L). We 
find a quite small, albeit significant, drift even in the largest lattices. Notice that the values 
for the ^/L ((74) crossing are monotonically increasing (decreasing) with L. Thus, unless 
something weird happens with the scaling corrections, Pc is bounded from above and below 
and one can readily extract pc = 0.311609(3). However, a more precise pc determination is 
possible using FSS techniques. We proceed as before to make a joint fit to Eq. (|^) for all 
the data in table ^, excluding those for L < Lmin- The results are presented in table |^. 

The value of w = 1.62(13) is remarkably larger than in most of the 3D systems (slightly 



L 


V 




.92 


94 


8 


0.8802(6) 


-0.01531(12) 


0.35395(11) 


0.72353(9) 


12 


0.8847(6) 


-0.03230(12) 


0.35395(11) 


0.72353(9) 


16 


0.8825(7) 


-0.03844(12) 


0.34854(10) 


0.72074(10) 


24 


0.8807(10) 


-0.04267(12) 


0.34601(10) 


0.71695(8) 


32 


0.8809(10) 


-0.04423(10) 


0.34559(10) 


0.71499(9) 


48 


0.8771(14) 


-0.04531(12) 


0.34603(11) 


0.71290(9) 


64 


0.8757(17) 


-0.04539(10) 


0.34638(10) 


0.71195(8) 


96 


0.8796(33) 


-0.04554(20) 


0.34672(24) 


0.71124(20) 


00 


0.8765(16)(2) 


-0.04602(27)(7) 


0.34675(26)(6) 


0.71052(21)(19) 


Table 5: 


Critical exponents for the site percolation obtained from pairs of type (L, 2L) 


obtained at Qf — 2. We also show the cumulants g4(i) and g2{L) at the same points. The 


last row 


correspond to an 


infinite volume extrapolation, showing the statistical error (first 


bar) and that coming from the uncertainty in u; 


(second bar). 




L 








Pc(.g4,i) 


8 




0.309761(7) 




0.313201(11) 


12 




0.311034(5) 




0.312454(7) 


16 




0.3113614(36) 




0.3120770(49) 


24 




0.3115337(20) 




0.3117950(26) 


32 




0.3115788(14) 




0.3117007(19) 


48 




0.3115992(9) 




0.3116390(12) 


64 




0.3116036(8) 




0.3116214(12) 


96 




0.3116063(9) 




0.3116122(11) 



Table 6: Crossing points for ^/L and 174 for the site percolation obtained from pairs of type 
(i,2i). 
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X'^/d.o.f. 




Vc 


24 


10.5/6 


1.57(2) 


0.3116092(5)(2) 


32 


0.65/4 


1.62(4) 


0.3116081(71(2) 


48 


0.07/2 


1.64(13) 


0.3116075(11) (2) 



Table 7: Infinite volume extrapolation of the crossing points for 54 and S,/ L, including data 
from L > Linin in the site percolation, in order to obtain to and Pc- 
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Figure 2: Qdp^ — 2^+" , Qx = 2^ '', 34 and 52 for pairs (L, 2L) as functions of X ".In the 
plot we use u = 1.62. 

below 1), but agrees with the prediction of e-expansions, and lies between the 1^20 = 2 of 
2D site percolation jl^] and the four dimensional value [0 uj4Y) = 1.13(10). Moreover, our 
value for cu is in agreement with other MC determinations (see table |^). 

In Fig. ^we display the quotients of x and dpS,, and the values of (74 and 32 as functions of 
L^^ for UJ — 1.62, measured at the points where = 2. The linear behaviour is much less 
clear than in the Ising case. One should not be surprised by this fact since for such a large w 
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Rcf. 


Pc 


cr 


T 


LO 




1 
1 


0.311604(6) 
0.311600(5) 


0.445(10) 


2.188(2) 
2.189(2) 
2.186(2) 


1.61(5) 

1.77(13) 

1.4 


This work 


0.3116081(7)(2) 


0.4522(8)(1) 


2.18906(6)(2) 


1.62(13) 



Table 8: Previous Monte Carlo determination of critical quantities for the 3D site percola- 
tion. 

it is very unlikely a clear separation between the leading corrections-to-scaling (as L~'^) and 
the sub-leading ones. One should be specially worried with the analytical corrections that 
for most operators go as L~^^^ « i^^. A parameterization of the sub-leading corrections is 
far from the present MC capacities. 

Fortunately, w is large enough to make the extrapolation almost unnecessary. For v, we 
do not find significant deviations for L > 48 and one could be tempted of simply averaging, 
obtaining ^'moan = 0.8768(10). However, we find no reason to consider as vanishing the 
coefficient of i^'^, and this assumption underestimates the errors. We find a non zero value 
of Aq^ in the fits to Eq. (^ for 77 and the cumulants. In the last raw of table || we present 
the results of these fits as well as the corresponding statistical errors, the second error bars 
corresponding to the uncertainty in uj. This w-error allows to quantify the possible shift 
that could be expected if the dominant corrections-to-scaling were the analytical ones, as 
the behaviour is basically linear with uj. One simply has to add 2.5 times the lo induced error 
to the central value for the extrapolation (the sign would be positive in the four cases). For v, 
rj and g2^ one can conclude that the systematic errors are hardly greater than the statistical 
one. For (74 the former could be twice the latter. Our final results can be contrasted with 
other MC estimates in table |^. 

6 Conclusions 

We have found that when measuring critical exponents and other universal quantities with 
high precision (below the 0.1%) with finite-size scaling techniques, a proper consideration 
of the corrections-to-scaling is mandatory. 

We have studied two simple three dimensional models. The Ising model shows corrections 
that can be parameterized with the leading corrections-to-scaling term. It is possible to 
obtain a very safe infinite volume extrapolation that can be as far as 10 standard deviations 
from the largest lattice's value. 

In the site percolation, the behaviour is completely different. The leading corrections- 
to-scaling cannot be easily isolated from the higher order ones, since the first irrelevant 
exponent is very large. However, its largeness makes the results on the largest lattices very 
near the infinite volume limit, and the difficulties of the extrapolation are not overwhelming. 
We have also measured with high precision the values of two universal cumulants (174, (72)- 
The non-vanishing value of the latter shows that the susceptibility is not a self-averaging 
quantity. 
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